Charge- versus spin driven stripe order: the role of transversal spin fluctuations. 
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The separation of the charge- and spin ordering temperatures of the stripe phase in cuprate 
superconductors has been used to argue that the striped phase is charge driven. Scaling analysis 
of a non-linear sigma model shows that the effect of spatial anisotropy on the transversal spin 
fluctuations is much more drastic at finite temperatures than at zero temperature. These results 
suggest that the spin fluctuations prohibit the spin system to condense at the charge ordering 
temperature, despite a possible dominance of charge-spin coupling in the longitudinal channel. 
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The observation of a novel type of electronic order in 
cuprate superconductors and other doped antiferromag- 
nets has attracted considerable attention recently. In this 
stripe phase, the carriers are confined to lines which are 
at the same time Ising domain walls in the Neel back- 
ground j^]. Substantial evidence exists that dynamical 
stripe correlations persist in the normal- and supercon- 
ducting states of the cuprates |2j . 

A further characterization of the fluctuation modes of 
the stripe phase is needed. In this regard, the finite tem- 
perature evolution of the static stripe phase might offer 
a clue. Both in cuprates [|l| and in nickelates the 
charge orders at a higher temperature than the spin, and 
both transitions appear to be of second order. Zachar, 
Emery and Kivelson |^ argue on basis of a Landau free 
energy that the stripe instability is charge driven: if the 
coupling between the charge- and longitudinal spin mode 
would dominate, charge and spin would order simultane- 
ously in a first order transition. This is a mean-field anal- 
ysis, and fiuctuations can change the picture drastically. 
For instance, at length scales larger than the interstripe 
distance the spin system remaining after the charge has 
ordered is just a quantum Heisenberg antiferromagnet in 
2-1-1 dimensions which cannot order at finite tempera- 
tures according to the Mermin- Wagner theorem. Zachar 
et al argue that the orientational ('transversal') fluctu- 
ations of the spin system can be neglected at the tem- 
peratures of interest, because it appears that the spin 
system left behind after the charge has ordered is not 
radically different from the antiferromagnet in the half- 
filled cuprates, exhibiting a Neel temperature of order 
300 K, an order of magnitude larger than that in the 
stripe phase. 

An important constraint is that the T = staggered 
magnetization in the stripe phase appears to be compa- 
rable to that at half filling j^. If the transversal fluc- 
tuations are responsible for the charge- and spin tran- 
sitions, it has to be demonstrated that the additional 
thermal fluctuations due to the presence of stripes have 
a much greater effect on the Neel state than the T = 



quantum fluctuations. To investigate this, we consider 
the simplest possible source of stripe induced spin disor- 
der. Following Castro Neto and Hone (CH) ^ we assume 
that the exchange coupling between spins separated by a 
charge stripe is weaker than the interdomain exchange, 
so that the collective spin fluctuations are described by 
a spatially anisotropic 0(3) quantum non-linear sigma 
(AQNLS) model. From our scaling analysis we find that 
a moderate anisotropy (a factor of ~ 4 difference in spin 
wave velocities) can explain a reduction of the Neel tem- 
perature by an order of magnitude, while the T = 
staggered magnetization is only reduced by a factor of 
two from its isotropic value. The reason can be in- 
ferred from the cross-over diagram (Fig. 1). As function 
of increasing anisotropy, the T = transition between 
the renormalized classical (RC) and quantum disordered 
(QD) states scales to smaller coupling constant, but the 
dimensionless temperature associated with the crossover 
renormalized classical- to quantum critical scales down 
much faster. Alternatively, we find that the behavior 
found by Chakravarty, Nelson and Ilalperin(CHN) for 
the correlation length in the RC regime of the isotropic 
model can be directly generalized to the anisotropic case: 
the expression for the classical anisotropic model remains 
valid when the bare stiffness is replaced by the renormal- 
ized stiffness. It is suspected that this holds more gen- 
erally. If so, the strong disordering influence of temper- 
ature as compared to the quantum fluctuation might be 
generic: whatever the disordering influence of the stripes 
is, it exerts it in an effectively three dimensional classi- 
cal system at zero temperature and in a two dimensional 
system in the finite temperature renormalized classical 
regime. 

It is assumed that the Neel order parameter fluctua- 
tions in the charge ordered stripe phase are governed by 
an AQNLS model [| [|, 

S AQNLS = 7^ [ dr f d^x (a {d^nf + {dyhf 
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where the bare coupling constant go and the spin-wave 
velocity c are those of the isotropic system, while a 
parametrizes the anisotropy. In the classical limit, this 
describes spin waves with velocity Cy{a) = c^/(Y+a}/2 
and Cx{a) = y/acy{a) in the y- and x directions, respec- 
tively. The slab thickness in the imaginary time direction 
u is given by phcA, where A is the cut-off of our spher- 
ical Brillouin zone. This model is derived by taking the 
naive continuum limit of a Heisenberg model with ex- 
change couplings J and a J in the y- and x directions, 
respectively. 

The renormalization of this model has received some 
attention recently [^[0]. We adopt here a variation on 
the procedure as proposed by Affleck 0. The central 
observation is that this model contains two ultraviolet 
cut-offs. As a ramification of the anisotropy, the highest 
momentum states in the x-direction will have an energy 
£jmax -^jiipji ig factor y/a smaller than that of the high- 
est momentum states in the y direction. Therefore, the 
initial renormalization flow from E^'^^ down to E^^^ is 
governed by one dimensional fluctuations. At E™^^ the 
resulting model can be rescaled to become isotropic, al- 
beit with 'bare' parameters which are dressed up by the 
one dimensional high energy fluctuations. 

Keeping the full model Eq. (|l|), the one dimensional 
fluctuations are integrated out (using momentum-shell 
renormalization |^ ) by neglecting the dispersions in the x 
direction entirely. This causes the anisotropy parameter 
a to become a running variable as well, which is always 
relevant. When the renormalized a = 1, the model has 
become isotropic, albeit with renormalized bare coupling 
constants. 

Writing h = (tt , a) , where a is the component of h 
in the direction of ordering, we expand to one-loop or- 
der in 7?. Subsequently, we Fourier transform the 7r-fields 
according to 



where 



Tf{x, t) = 



IK-X — lUJnT 



(2) 
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where w„ = 27m/u are the Matsubara frequencies. The 
momenta k are rescaled with A to become dimensionless. 
Separating the flelds according to 



7f(A:, n) 



7f>(fc,7l) 



< |fc„l < 1 



7f< (fc, n) ; < \ky \ < e 
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where I is small, we integrate out the flelds 7r>, using a 
square Brillouin zone for convenience. Rescaling ky, 7r<, 
u, g and a, we find that the model scales to larger a 
(smaller anisotropy) . We obtain the following flow equa- 
tions 



a = aoe , 

d£ ^ a 

dl 
dt 
dl 



l + a 
t + tgl, 



9 + 9^1, 
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^ VTTT^ fi^^ coth(ly^V^fcgTT) 

4V27r2 " ^akl + 1 
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and where t is the dimensionless temperature, to = 
fceT/pg. From eq. (|) and (|), wc find for the slab 
thickness u = g/t 
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From Eq. (Q) it follows that a = 1 corresponds 
with ^ = ^1 = — In y^oq. The bare coupling constant 
(at T = 0) and bare slab thickness of the effective 
isotropic model follow by integrating Eq.'s (||j6|) down 
to h {gi = g{li),ti = t{li)), 



91 = go/ 
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(arsinh(ya^^)/ya^^- 



1 + Qfo 27r2 
ln(l + VTT^) - ln(V^(l + . (9) 

51 Ai = (goAo)Vao(l + ao)/2 (10) 

Except for these altered bare quantities, the isotropic 
model is analyzed in the standard way Q . 

Putting gi = gc — 47r and solving go, we find the crit- 
ical bare coupling for the anisotropic model 



1 + ao 



(arsinh(^/a^^)/Vc^o + 



ln(l -f yrr^) - ln(V^(l + V^)2))] . (11) 



Wc find this result to be the same within a couple of 
percents as the outcome of large N mean field theory 
while the difference originates in an inaccuracy in 
our calculation related to the switch from the square (at 
E > i^iJiax) to the spherical Brillouin zone of the effec- 
tively isotropic model. 

For ao — I, the one- loop cross-over lines between 
the QC and the RC/QD regime are given hy t ~ 
±27r(l — g/47r). Taking (51, ti) to lie on these lines and 
iterating the flow equations backwards, we obtain the 
cross-over diagram for the anisotropic model, shown in 
fig. |l|. Note that the anisotropy has a stronger effect on 
the i-dependence of the RC to QD line than on its g- 
dependence. This already indicates that the T = prop- 
erties will be less affected by the anisotropy than those 
at finite temperatures. 
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FIG. 1. Cross-over diagram for the anisotropic QNLS. The 
lines are for a = 1, 0.4, 0.1 and 0.025 from top to bottom. The 
end-points of the quantum-critical to quantum-disordered 
lines map onto = (87r,27r). Notice that when to be- 

comes larger than the crossover temperature from renormal- 
ized classical to quantum critical at go = one dimensional 
fluctuations are dominating for all values of go- 

The one-loop mapping to an isotropic QNLS provides 
a simple way of calculating the correlation length in the 
anisotropic model. Noting that the correlation length in 
the y-direction scales as ^ = ^oe~' under Eq. (||), it imme- 
diately follows that £,{go,to) = e'i^isotr.(5i, ^i)- Inserting 
the 1-loop expression for ^isotr. in the RC regime |^ and 
using Eq.'s (|9|,p^) (the use of the T = expression for gi 
Eq.(p|) is a good approximation if gi/ti ^ 1), 
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exp[V^p,(0)/fci3T]. (12) 
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where the renormalized T = stiffness is given by, 
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Eq.'s ( |l2yi3| ) is our central result. It shows that the cor- 
relation length in the renormalized classical regime has a 
twofold exponential dependence on the anisotropy, both 
originating in the high frequency one dimensional fluctu- 
ations. As already pointed out by CH ||^, the anisotropy 
causes gc to decrease (e.g.. Fig. 1), leading to a reduc- 
tion of ^ at a given temperature. However, we find an 
additional y/a in the exponent which has been overlooked 



before, and this is the specific way in which the greater ef- 
fect of the thermal fluctuations as we noted earlier shows 
up in the renormalized classical regime. In fact, it shows 
that the basic invention of CHN |6| is straightforwardly 
extended to the anisotropic case. The correlation length 
is given by the expression for the classical system, and 
quantum mechanics only enters in the form of a redefi- 
nition of the stiffness. However, for the classical corre- 
lation length expression one should use the one for the 
anisotropic classical model. Using the same procedure as 
for the quantum model, it is easy to demonstrate that 
the correlation length of the anisotropic classical 0(3) 
model in 2D behaves as ^ exp{y/aop^/kBT), and this 
explains the occurrence of the additional y/oo factor. 

The finiteness of the Neel temperature is caused by 
small intraplanar spin- anisotropics and interplanar cou- 
plings. Keimer et al |^ have shown that in La2Cu04 the 
former dominate, and these can be lumbed together in a 
single term aef / which plays the role of an effective stag- 
gered field. The Neel temperature can be estimated by 
comparing the thermal energy k^T^ to the energy-cost 
of flipping all spins in a region the size of the correlation 
length in the presence of the effective staggered field. 
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a Mo 
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Because it is not expected that stripes will influence 
the spin anisotropics strongly, we can use the estimate 
for Q!off as determined for the half-filled system: a^s = 
6.5 X 10"'' 1^. For our estimate of T/v, we will use 
spin- wave results for the renormalized stiffness, suscepti- 
bility and spin wave velocity [Q. For S — 1/2, they 
are tic = 0.5897v^Ja, X-l(O) = 0.514?iV(8 Ja^), and 
Ps = c^x_l(0). The bare coupling constant is obtained 
from (.9o/4vr) = 1/(1 + A-Kx^c/hh) [§, which yields 
(7o — 9.107 for Aa — 2^pK. We notice that the 1-loop 
result for the prefactor is not correct, but this factor is 
not very important as far as the reduction of the Neel 
temperature is concerned. 

Since our T = results coincide with those obtained by 
CH , we use their expression for the zero temperature 
staggered magnetization 1113] , 



M,(a) 



l-.9o/47r ' 
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and its anisotropy dependence is shown together with the 
results for the Neel temperature in the inset in fig. |[ To 
illustrate the effects of a different aoff in the stripe phase 
(e.g., Jj^ may be much reduced due to frustration) we 
have also plotted the results for a^sioL < 1) = 10acfr(a = 
1) (upper dashed line) and for aoff (o; < 1) = 0.1acfr(a — 
1) (lower dashed line). In Fig. ^T^v is plotted versus Mg. 
As expected, the dependence of Tjv on anisotropy is con- 
siderably stronger than that of Mg ■ A reduction of Ms by 
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a factor of 2 due to a spin- wave anisotropy ~ yja ^ 1/A 
order is accompanied by a suppression of T/v by roughly 
an order of magnitude. 
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FIG. 2. The Neel temperature versus the zero temper- 
ature staggered magnetization, with the anisotropy as im- 
plicit parameter. Both quantities are normahzed with re- 
spect to their value in the isotropic system. The upper/lower 
dashed line gives Tjv for aoff (Neel stabilizing field) a factor 
10 larger/smaller than in the isotropic system. Inset: Tn and 
Ms as a function of the anisotropy parameter a. 

In the above we relate different experimentally acces- 
sible quantities (spatial- and spin anisotropics, Neel tem- 
perature, T = staggered order, correlation length) and 
further ]l3| , [r^ experimentation is needed to unambigu- 
ously demonstrate that spatial anisotropy is the cause 
of the low spin ordering temperature. If the fluctuation 
behavior in the RC regime is indeed as general as sug- 
gested by the present analysis, other sources of stripe in- 
duced spin disorder could have similar consequences. For 
instance, local charge deficiencies in the stripes caused 
by quenched disorder would give rise to unscreened (by 
charge) pieces of domain walls. Such stripe defects are 
like the dipolar defects discussed by Aharony et al, p^ ] 
and their frustrating effect is expected to be dispropor- 
tionally stronger at finite temperature than at zero tem- 
perature. 

Above all, the present analysis shows that a Landau 
mean-field analysis falls short as a description for the 
thermodynamic behavior of the stripe phase because of 
the importance of fluctuations. Stronger, thermodynam- 
ics does not offer an unambiguous guidance regarding 
the microscopy (frustrated phase separation |l6[ versus 
'holon' type mechanisms |]l^). Here we have focussed on 
the transversal spin fluctuations, and given that there is 



ample evidence for a pronounced slowing down of the spin 
dynamics at the charge ordering temperature, these un- 
doubtedly play an important role. It is noted that recent 
results point at a similarly important role of fluctuations 
in the charge sector [Q. 
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